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Abstract 

We use an appropriate combination of moments of finite energy sum 
rules in QCD in order to compute the _Bg-meson decay constants fs 
and /b, -We perform the calculation using a two-loop computation of 
the imaginary part of the pseudoscalar two point function in terms 
of the running bottom quark mass. The results are stable with the 
so called QCD duality threshold and they are in agreement with the 
estimates obtained from Borel transform QCD sum rules and lattice 
computations. 
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1 Introduction 



Since the pioneering work of Shifman, Vainshtein and Zacharov T, Laplace 
sum rules have been successfully applied to calculate all sort of parameters of 
the hadronic spectrum. The main advantage of this type of sum rules is that 
a Borel transform applied to the correlation function enhances the contribution 
of the low lying resonances of the hadronic spectrum, which properties are to 
be determined. On the other hand, it reinforces the convergence of the QCD 
asymptotic calculation in the high energy domain. The price to pay with this 
method is the appearance of the so-called Borel parameter which has to be fixed 
by stability arguments. Other particular sum rules based on Hilbert transform, 
inverse moments.., have also been used to suitably deal with other particular 
problems. 

In this note we introduce a method based on positive moments of QCD finite 
energy sum rules. Traditionally, this type of sum rules have the disadvantage 
of reducing the contribution of the low energy part of the hadronic spectrum, 
whereas they enhance the QCD high energy region. Although they are easy to 
handle, one needs to fix a QCD duality threshold, where theoretical calculations 
are accurate enough and at the same time the low energy region admits a suitable 
hadronic parametrization. Nevertheless, it is not always easy to fix the value of 
this duality threshold in the sense that the results have to be independent of this 
value. The method we propose is to combine different moments of finite energy 
sum rules in order to get a polynomial weight for the correlation function such 
that its contribution becomes negligible when integrated from the continuous 
physical threshold until the QCD duality threshold. In this way, the low lying 
resonance region is enhanced in the hadronic spectrum and, on the other hand, 
the role of the QCD duality threshold becomes less relevant on the final results. 
We will be more explicit on this in the next section. 

Here we use our method to evaluate the decay constants of the lightest 
pseudoscalar bottom mesons (/s and /s,), which parametrize the i?g-meson 
matrix elements of the pseudoscalar current: 

<n\{Mb + mg)qi-f5b)iO)\B, > = /b,M|^. 

These decay constants have received recently a lot of attention since they 
enter in hadronic matrix elements oi B — B mixing, and its accurate evaluation 
would facilitate a better determination of the Bb mixing factor from recent 
experiments carried out in the B-factories. They also appear in the leptonic B 
decay widths and its knowledge could provide a good determination of the \ Vqb\ 
matrix element in future experiments. Calculations of these decay constants 
have been performed since the eighties, with results in the range of fs = 160 — 
210 MeV and fsJlB = 1-09 - 1.22 from Borel transform techniques H H El . 
Computations in lattice QCD give also results in a wide range: fs = 161 — 
218 MeV and JeJIb = l-H - 1-16 [3 El (for a review and a collection of the 
results, see [7]). As we see, in these range of values there is still room for 
improvement. 



2 



The plan of this note is the following: in the next section we briefly review 
the theoretical method proposed, in the third section we discuss the theoretical 
and experimental inputs used in the calculation and in the fourth one we present 
our results for the decay constants with a discussion of the errors. We finish the 
paper giving the conclusions. 



2 The method 

The two point function relevant to our problem is: 

n(s = q^) = ij dx e*"^ < n\T{j5{x)j5m \n >, 

where < f2 1 is the physical vacuum and the current j5{x) is the divergence of 
the axial- vector current: 

is (a;) = {Mq + ruq) : q{x) 175 Q{x) : 

Mq is the mass of the heavy quark Q{x) which will the bottom quark in our 
case, whereas stands for the light quark mass, up, down or strange. In order 
to write down the sum rules relevant to our calculation, we apply Cauchy's 
theorem to the two point correlation function n(s), weighted with a polynomial 
P{s) as indicated in the following: 

^^£s'P{s)U{s) ds = 0, (1) 

(the power s% with i > 0, is introduced here for convenience, as it will become 
apparent in the following). 

The integration path T is extended along a circle of radius \s\ = sq, and 
along both sides of the physical cut starting at the physical threshold Sph., i- 
e. running in the interval s € [sph., so]- Neither the polynomial P{s) nor the 
power of the integration variable change the analytical properties of H{s), so 
that we obtain the following sum rule: 



-/ 



" s' P(s) Im n(s) ds = - ^ (f s' P{s) n(s) ds (2) 



On the left hand side of this equation, we enter the experimental information 

of Imn(s), starting from the physical threshold Sph. till the integration radius 
So, whereas, on the right hand side, we consider the asymptotic QCD theoretical 
calculations to be plugged into the integration contour of radius so- Therefore, 
this radius makes the compromise where the contribution of the hadronic spec- 
trum can be approximated by the QCD calculation. This is our QCD duality 
threshold that we refer to in the previous section. 
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The asymptotic expansion of QCD (11'^*-'^ (s)) can be split in two parts, 
including the perturbative and non-perturbative terms, as follows: 

nQCD(g) = np"*(s) + n"°"p"*(s), (3) 

At this stage, we consider that nP"* (s) is an analytic function of s, with 
a real cut starting at sqcd — (Mf, + m^) , therefore, we can use again the 
Cauchy's theorem to convert the integration of nP°''* (s) along the circle \s\ = Sq 
into an integration of the corresponding absorptive part along the QCD cut. 



/P(s)lmn(s) ds = 

'P{s)lmW'' {s) ds-^ I s'P(s)n"°"P"* (s) ds 
27rz J|s|=so 



(4) 



On the right hand side of equation we take for the perturbative spectral 
function the exact two-loop QCD calculation [J 
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where Mf, = M{, (/i) is the running mass of the b quark in the MS scheme. It is 
known that the expansion in terms of the running mass converges much faster, 
in a wide range of the renormalization scale, than the one in terms of the bottom 
pole mass, as notice, for instance, in (4j. 

For the non perturbative part n"°"P'^''* (s) we take the contribution com- 
ing from the vacuum expectation values of non perturbative operators up to 
dimension six E] ■ 



jj„onport.(^g) = 



Mpoic <qq> < —G^ > 



- -M: 



pole 



(5- 



(^-Mp'ole)' 



< qaGq > - 
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(^-A^p'ole)' (.'^-M^^le)^ (S-M2„,e)^ 



(6) 



The contour integration of the non perturbative part is easily done by means 
of the residues theorem. 

Again, for calculational purposes, we consider in this non perturbative ex- 
pansion the relation between the pole mass Mpoie and the running mass Mb in 
the appropriate order of the couphng constant ^ El E] • 

Finally, on the left hand side of equation we parametrize the absorptive 
part of the two point correlation function by means of a narrow width approx- 
imation of the lightest Bq resonance plus the hadronic continuum of the bq 
channel starting at Scont, above the resonance region: 

- Imn(s) = fl Sis - A/|J + ilmn-"'-0(.s - w.) (7) 

where Mb and Jb are respectively the mass and the decay constants of the 
lowest lying pseudoscalar meson Bq. 

Looking back to equation (O and taking into account all the theoretical 
parameters as well as the mass of the S^-meson as our inputs in the calculation, 
we see that the decay constant can be computed as far as we had a good control 
of the hadronic continuum contribution of the experimental side. 

Since this is not the case, to cope with this problem we make an appropriate 
choice of the polynomial {P{s)) in that equation. We take: 

P{s) =ao + ais + oas^ + a^s^ + ... + a^s"^ , (8) 

such that its coefficients are fixed by imposing a normalization condition at 
threshold P(sph. = M^^) = 1, and requiring that should vanish in the range 
[scont.jSo] in a least square sense, i.e.. 

So 

s'= P{s) ds = for = 0, . . . n (9) 

cont. 

These conditions exactly cancel the continuum contribution as far as Im n'^""*- 
can be well approximated by an n degree polynomial. On the other hand, 
by virtue of the normalization condition, it will enhance the role of the Bq 
resonance. Notice however that by increasing the degree of the polynomial, 
P(s), we will require the knowledge of further terms in the non perturbative 
series, which are unknown. Therefore, a compromise criterion on the choice of 
the polynomial degree has to be taken. 

To check the consistency of the method in this work, we have considered a 
second and third degree polynomials and the results are fully compatible within 
the range of the errors introduced by the inputs of the calculation. We also 
have checked explicitly that a smooth continuum contribution has no influence 
in the result. This procedure was previously used in the calculation of the 
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charm mass from the cc experimental data. The continuum data from the 
BES II coUaboration |13| had no influence when an appropriate polynomial was 
included Employing the same technique, very accurate prediction of the 

bottom quark mass was also obtained using the experimental information of the 
upsilon system [TH] , 

After these considerations we proceed with the analytical calculation of the 
decay constant /b, , neglecting the contribution that comes from Imll'^""*-. For 
different powers of we obtain different sum rules which, in principle, should 
give the same result for the decay constant Jb^ ■ Notice however that, by increas- 
ing the value of i we get larger contribution from both, the large s region of the 
spectral function and higher orders in the non perturbative QCD expansion^. 
Since these pieces have the main experimental and theoretical uncertainties, we 
will consider the calculation coming from the sum rules corresponding to the 
powers with i = and i — 1 which, in any case it is enough for our purposes. 

M^B, fl^ = - r P(s)ImnP-' (5) ds -Res {P(s)n"°"P"* (s), s = M^} 

(10) 

for i = ("first sum rule") and 

M% fl = - f " sP(s)ImnP"' (s)ds-Res|sP(s)n"°"P°^* (s), s = M?] 

(11) 

for 1 = 1 ("second sum rule"). Let us emphasize that in these two sum rules 
there are two unknowns, the decay constant Jb and the QCD duality threshold 
So, the last one appears in the upper integration limit and also in the coefficients 
of the polynomial P(s). Therefore, we can use both sum rules to determine /s^ 
as well as sq- To employ a couple of sum rules in order to fix the QCD duality 
threshold (sq) is a usual procedure [21 El, but, as we will see, with our method 
we have the additional advantage that the value obtained for sq is very stable, 
in the sense that any change of this value would not affect appreciably the result 
of the decay constant fs^ ■ 

3 Results 

In the calculation of Jb we take = in the factor {Mb + ruq)'^ of the 
correlation function, and in the low integration limit of equation Q , whereas for 
fss we keep the strange quark mass different from zero although its contribution 
turns out to be negligible. 

The experimental inputs are as follows. For q being the light quarks q = u, d 
we have the physical threshold Sph., at the squared mass of the lowest lying 
resonance in the bq channel : 

Sph. = Ml = 5.279^ GeV^ (12) 

■^The results obtained with higher moments are compatible with the ones found here. 
However they are less stable with the duality parameter sq and, hence, less accurate. 
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The continuum threshold Scont. is taken at the Btttt intermediate state, in an 
/ = ^ s-wave, i. e. 

Scont. = [Mb + 2m^f = 30.90GeV^ . (13) 

For the strange s quark we take: 

Sph. = Ml^= 5.369^ GeV^ 
Scont. = (Mb, + 2to,)' = 31.92 GeV^ (14) 

In the theoretical side of the sum rule the inputs we take are as follows. The 
strong coupling constant is taken at the scale of the electroweak Z boson mass 

ini 

a^Afz) = 0.118 ±0.003 (15) 

and run down to the computation scale using the four loop formulas of reference 
|16| . The values, with the corresponding error bars, for the quark and gluon 
condensates (see for example 0]) and the mass of the strange quark are: 

(gg)(2GeV) = (-267 ± 17 MeV)^ 
i—GG) = 0.024 ± 0.012 GeV, 

TT 

{qcrGq) = ml{qq), (with = 0.8 ± 0.2 GeV) 

(ss) = (0.8±0.3)(gg) 

TO, (2 GeV) = 120±50MeV (16) 

We also need to fix the renormalization scale /i in the equations |(SJ) and For 
definiteness we take /i = Mft(Mfc), although in order to find the stability of the 
results under the renormalization scale we will study its variation for values of 
this scale in a range appropriate for our discussion, namely, (Aq^^ ^ < sg). 

Finally for the bottom quark, the value Mi, « 4.20 GeV is nowadays generally 
accepted. We take the result of 15 which has also been obtained with the sum 
rule method described here: Mb(Mb) = 4.19 ± 0.05 GeV. 

Now, we proceed in the way described before. Firstly, we compute/s as a 
fimction of sq with the two different sum rules for i — 0,1. Then, we fix so 
at the value where both sum rules give the same result for fs- With this usual 
procedure, and for a second degree polynomial, we find so = 48.5 MeV^ and 
fs = 185 MeV. 

From the theoretical inputs quoted above the main source of errors comes 
from the bottom mass which, in the range given above, produces a variation 
in the decay constant of ^ 11 MeV. Other source of errors are the quark con- 
densates"^ which affect the result by ±4 MeV. Since the decay constants are 
physical observables, the results should be independent of the renormalization 

^As for the contribution of the condensates, the only relevant comes from the lowest di- 
mension one which gives an 8% of the total result. The higher dimension terms considered 
here doesn't give a sizeble contribution giving a hint on the convergence of the condensate 
series in this approach. 
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Figure 1: Decay constant Jb as a function of the integration radius sq for 
Mi,{Mb) — 4.19 GeV. With a second degree polynomial in the sum rule the 
dashed line represents the case i — 1 and the solid line the case z = 0. 



scale /i. However, here we have used an approximation in the asymptotic spec- 
tral function of QCD, taken only the two- loop order (|5I6|I . Therefore a residual 
dependence on the renormalization scale is expected. In order to quantify the 
uncertainty of fixing the scale at the bottom mass, we vary the scale in the range 
H e [3,6] GeV, introducing an uncertainty of ±18MeV in the result. This de- 
pendence in the renormalization scale is expected to lower down if higher orders 
in the coupling constant in (|5I6|I are taken into account. 

Adding quadratically all these errors, we finally quote the following result 
for the decay constant of the light meson B: 

fs = 185 ± 22 MeV. (17) 

Notice in figure 1 that with i = 1 (dashed line) and a second degree poly- 
nomial there is a stable value (minimum), /s — 183 MeV, at sq — 44GeV^. 
But for i — (solid line) the result Jb — 180 MeV, is an stable value (infiexion 
point) which is practically constant around sq = 60 GeV^ However, to compare 
with other results from QCD sum rules, we take the crossing point of figure 1 
as our final result. The stable results of both curves are completely compatible 
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within error bars in such a way that we could have considered the value for /s 
in the wide range of the stability region of sq. 

Proceeding in the same fashion for the Bs meson, with the only change given 
by the nonzero mass of the light s quark, we find the decay constant /b, ■ In 
this case, as can be appreciated in figure 2, the Sq value where the two sum rules 
give the same result for /b^ is sq = 49.6 GeV^. 




Figure 2: Decay constant fs^ as a function of the integration radius sq for 
Mi,{Mb) — 4.19 GcV. With a second degree polynomial in the sum rule the 
dashed line represents the case i = 1 and the solid line the case z = 0. 

The result in the intersection point is: 

Jb, = 202 ± 24 MeV. (18) 

where a similar analysis of errors has been considered. The only new ingredient 
is the uncertainty coming from the strange quark mass, which contribution turns 
out to be negligible. 

It is of special interest the ratio of the decay constants fs^ and fs, which 
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should be 1 in the chhal hmit. We find 

^ = 1.09 ± 0.01 (19) 

JB 

We are free to remark that in the calculation of this ratio, the uncertainties 
of the theoretical parameters are correlated, this is why the final error becomes 
very small. 



4 Conclusions 

In this note we have computed the decay constant of _Bg-mesons for q either the 
strange s or the u or d massless quarks. We have used a suitable combination of 
moments of QCD finite energy sum rules in order to minimize the shortcomings 
of the available experimental data. On the theoretical side of the sum rule, 
we have used the pseudoscalar two point function calculated up to two-loop in 
perturbative QCD and up to dimension six condensates in the non perturbative 
QCD expansion. Instead of the commonly adopted pole mass of the bottom 
quark, we use the running mass to get a good convergence of the perturbative 
series. We have a good control of the results against the duality threshold Sq, 
which turn to be very stable. 

The results found, taking the running mass of the bottom quark Mi,{Mi,) = 
4.19 ± 0.05 GeV are given in equations (|17I18|) and we collect here for conve- 
nience: 

fs = 185 ± 22 MeV, /b, = 202 ± 24 MeV. 

where the error bars come from the uncertainty in the theoretical parameters 
H16(l as well as the residual dependence on the renormalization scale. We notice 
that the results are very sensitive to the value of the running mass, giving most 
of the theoretical uncertainty. On the other hand they turn out to be quite 
stable against the variations of the other input parameters, in particular the 
integration radius Sq. 

In this treatment we could not include the three-loop corrections in the two 
point correlation function, since the complete analytical QCD expression along 
the cut is not known. Despite of this lack of information, one can interpolate the 
three- loop low energy QCD calculations with the high energy ones 0]^]. The 
results do not diff'er much from the exact two-loop calculation. Another possi- 
bility, is to include the three-loop high energy expansion along the integration 
circle \s\ — so which is known to certain accuracy. We have also performed this 
calculation somewhere else, having results very close to the ones found here. 

Within the error bars, our results agree with other results in the literature, 
obtained with different sum rule methods and with lattice computations |21 01 
El mini- However, we advocate for values of the decay constants in the lower 
band. 
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